Emergent p-wave Kondo Coupling in Multi-Orbital Bands 
with Mirror Symmetry Breaking 
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We examine Kondo effect in the periodic Anderson model for which the conduction band is of 
multi-orbital character and subject to mirror symmetry breaking field imposed externally. Taking 
p-orbital-based toy model for analysis, we find the Kondo pairing symmetry of p-wave character 
emerges self-consistently over some regions of parameter space and filling factor even though only 
the on-site Kondo hybridization is assumed in the microscopic Hamiltonian. The band structure 
in the Kondo-hybridized phase becomes nematic, with only two-fold symmetry, due to the p-wave 
Kondo coupling. The reduced symmetry should be readily observable in spectroscopic or transport 
measurements for heavy fermion system in a multilayer environment such as successfully grown 
recently. 
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Introduction.- State-of-the-art technique to grow heavy 
fermion compounds with atomic thickness precision has 
been achieved in recent years [1, 2]. Heavy fermion phase 
and superconductivity were identified in CeCoIn 5 layer of 
less than one nanometer thickness grown in this manner, 
prompting speculations of exotic new phases in atomi- 
cally thin heterostructure [2]. Apart from the novelty of 
reduced dimensionality, the alternate stacking of heavy 
fermion layer with normal metal YbCoIns layer realized 
in the growth process is likely to give rise to new physical 
phenomena by virtue of the mirror symmetry breaking 
(MSB) at the interface. Aspects of non-centrosymmetric 
superconductivity have been exhaustively studied since 
its original discovery in the heavy-fermion compound 
CePtaSi [3]. The heavy fermion multilayer structure, 
on the other hand, differs from the conventional non- 
centrosymmetric material in that the violation of sym- 
metry is macroscopically imposed (not microscopically at 
the level of crystal structure) and, with recently available 
techniques such as electrolyte gating, controllable. Due 
to the thin layer thickness, such MSB is likely to infect 
the entire heavy fermion band structure and therefore 
the nature of emerging Kondo and/or superconducting 
pairing as well. Some consequences of the MSB in the 
heavy fermion multilayer have been studied theoretically 
by Maruyama et al quite recently [4], where the main 
focus was on the novel superconducting pairing symme- 
tries and other thermodynamic properties induced by the 
phenomenological Rashba term. 

It has been customary in the past to regard the Rashba 
interaction as a relativistic effect manifesting itself in a 
mirror symmetry-broken environment. Several recent pa- 
pers [5-7], however, have shown that multi-orbital degen- 
eracy, in addition to MSB, is an essential pre-requisite 
for the manifestation of strong Rashba coupling in the 
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FIG. 1: (Color online) Schematic picture of Kondo lattice 
model subject to mirror symmetry breaking at the surface or 
interface. Hopping processes between, for instance, p z and 
p x , y otherwise forbidden in the bulk become possible due to 
the symmetry breaking field as shown in the upper right car- 
toon. 



electron dynamics. In this new picture Rashba interac- 
tion becomes an electrostatic effect with its magnitude 
dictated by the degree of orbital asymmetry about the 
mirror plane (see Fig. 1) [5-7]. 

In fact Kondo lattice behavior takes place in an 
orbital-rich environment with the extended p- or d- 
orbital states forming the conduction bands while the 
localized moments arise from /-orbit als [8, 9]. Using 
our recent experience with the Rashba phenomena in 
non-Kondo materials as a guide, we present a faithful 
modeling of Kondo lattice behavior where multi-orbital 
degrees of freedom in a mirror symmetry-broken en- 
vironment are explicitly embedded in the conduction 
band Hamiltonian. The use of Rashba interaction as a 
purely phenomenological addition, as is often done in 
the current literature, is avoided. Our approach differs 
as well from the large- TV generalization of the one-band 
Kondo lattice model, where N identical copies with the 
same orbital symmetry are adopted as a way to justify 
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the mean- field approximation [9]. In reality the MSB 
field - usually the electric field acting perpendicular to 
the interface or the surface - influences orbitals which 
are extended along the field direction more than those 
which extend in the orthogonal direction. Conventional 
one-band or large-iV approach assuming featureless 
orbitals will not capture this effect properly. 



The model- Choosing the degenerate p-orbital states 
hopping on a square lattice as a toy example, we write 
down the Hamiltonian [6, 7] 



HuSB = ( e kPAka^Aka + SjiSxP^PX^a + h.C.) 
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FIG. 2: (Color online) Phase diagrams with (a) a = (no 
SOC), n c = 2 and (b) a = 0.5 (SOC), n c = 1.5. Phase I 
denotes all three p-orbitals hybridized (h%, h v a , h% ^ 0). Phase 
II represents only p y orbitals hybridized (h% ^ 0, /i^ 7^ 0, 
/ij = 0). Phase III is trivial with no Kondo hybridization 
(h& — 0). In (a), as shown by a red vertical arrow, we may 
encounter two distinct Kondo temperatures labeled by X£ 
and T^ y . Parameters used to obtain both phase diagrams are 
(Vi, V2, Jk) = (2, 0.65, 3). Fermi surface structures at the four 
points indicated by stars are shown in Fig. 3. 



for the dispersive p-orbital (p\k a ) bands and a single pair 
of localized Kramers doublet (fwa) with on-site energy 
/if. The orbital index A G {x,y,z} spans the three p- 
orbital states of the conduction band and a = +(— ) is 
the spin up(down). One can show, employing the two 
Slater-Koster parameters Vi, V2 for a- and 7r-bonding 
of p-orbitals, that e£ = -2(Vi + V 2 )c\ + 2V 2 (c x + c y ) 



with c x ( y ) = cos/c. 



x(y)i s x{y) 



— sin/c. 



and c z = s z = 



for the assumed two-dimensional bands. Mirror symme- 
try breaking is reflected in the second term of Eq. (1), 
proportional to the MSB parmaeter 7 [7] . They give the 
hybridization of p^-orbital with p x - and p^-orbitals; these 
are processes normally forbidden by symmetry if a mirror 
plane existed. This kind of symmetry breaking is realized 
naturally on the surface of a bulk material or by applying 
vertical electric field to the planar system as illustrated 
in Fig. 1. 

Secondly the Kondo hybridization of local and itin- 
erant spins follows from the Heisenberg-type exchange 
interaction Hk = Jk 52i a • Sif, with spin- 1/2 oper- 
ators made out of £>A-orbital (Si\) and /-orbital (S^/), 
respectively. Finally i^soc = ( a /2) J2i • en involving 
the p-orbital angular momentum operator (L = 1) 
and the spin- 1/2 operator (1/2) 0^ at each site i gives the 
spin-orbit coupling (SOC) within the p-orbital bands. In 
typical heavy fermion matter such as CeCoIns the heavy 
element In playing a main role in the conduction band 
structure will also carry substantial amount of spin-orbit 
energy a. 

The total Hamiltonian H = i^MSB + ^soc + ^K can be 
handled within the framework of renormalized mean- field 
theory [10-12] which decomposes the Kondo exchange 
interaction as (a = —a) 



+ yE ? p jU. + J -f £ r^Uw (2) 

i<7 iXcr 

In addition to the Kondo hybridization parameters 
b x aa , = {flpxi*') and h x = b x a + 6^/2, we also 
allow same-orbital spin-flip processes encoded by 

T l = T,\(P\iaPXi(t) and T t = (fiafi*)- The Spin-flip 

parameters r£, and can be nonzero due to the 
spin-orbit interaction. We restrict ourselves however to 
non-magnetic, singlet Kondo coupling in the mean- field 
scheme and avoid the possibility of Kondo pairing be- 
tween remote sites [10], which is another known avenue 
to obtain non-s-wave Kondo pairing. The fact that our 
model can exhibit the p-wave Kondo phase despite the 
on-site-only Kondo hybridization Hamiltonian Hk em- 
phasizes the emergent nature of unconventional pairing. 
As demonstrated in Refs. [5-7], Rashba interaction is 
another emergent property of the band that can be 
derived from the underlying microscopic Hamiltonian 
consisting only of Hmsb and i^soc- 

Phase diagram - All the hybridization parameters can 
be determined self-consistently, along with /if to ensure 

(52a fiafiv) = 1' ano - ^he overall chemical potential /i, 
for various input band structure parameters Vi,V2,7, 
electron filling factor n c = (52\aP\iaPM<r)i a ^ various 
temperatures T. Generic phase diagrams resulting from 
our model Hamiltonian are plotted in Fig. 2 in the space 
of temperature T and MSB parameter 7 with and with- 
out the SOC. We observe two distinct Kondo-coupled 
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FIG. 3: (Color online) Band dispersions at the Fermi level 
for the various phases marked by red stars (★) in Fig. 2. (a), 
(b) and (c) are the Fermi surfaces of the phases I, II and 
III in Fig. 2(a) respectively, and (d) corresponds to the star 
in Fig. 2(b). The reduced symmetry in the Kondo-coupled 
phases (a), (b), and (d) due to the p-wave Kondo pairing are 
apparent against the (k x , k y ) axes in red dashed lines. 



phases labeled I in Fig. 2, where all three p-orbitals 
hybridize with /-orbital ft,*/, ft* ^ 0), and II where 
only the in-plane p x ,p y orbitals hybridize (ft* = 0). 
The high-temperature non-Kondo phase with = 
is labeled III. Both transitions, I — ^11 and II— >III, are 
second-order. Phase I becomes unstable and gives way 
to phase II as shown in Fig. 2(b) as the strength of SOC 
is increased. Although the addition of SOC complicates 
the microscopic Hamiltonian, fortunately, many of the 
hybridization parameters were found to be vanishing 
Tl = Tl = b z aa = b x a - a = b% = even with finite SOC. 
Features of the phase diagram are dependent on the 
filling factor n c of p-electrons. Phase I is observed within 
a small range of n c and, as indicated in Fig. 2(a), may 
even exhibit some re-entrant behavior with increasing 
7. While the details regarding the extent of each Kondo 
phase or the Kondo temperature are sensitive to the 
occupation n c , the effect of SOC on them was found to 
be much less significant by comparison. 

Effective Hamiltonian for p-wave Kondo pairing.- Most 
interestingly, the two Kondo-coupled phases I and II ex- 
hibit the band dispersion with the loss of symmetry un- 
der the reflection (k Xl k y ) — > (—k Xl k y ) or (k Xl k y ) — >• 
(k x ,—k y ), as well as the 90° rotation, as shown in Fig. 
3. Reduction in the band structure symmetry only oc- 
curs in the Kondo phase implying that it can only be 
the consequence of p-wave Kondo coupling. Below we 
show rigorously how such unconventional pairing occurs 
for MSB-infected bands. 



With a proper choice of tight binding parameters and 
filling factors we can identify cases in the self-consistent 
band structure wherein the two lowest-energy spin-orbit 
coupled p-orbital bands and the two Kramers doublet f- 
electron bands occur in the vicinity of the Fermi energy 
while all others are further removed from it. In such cases 
one can describe the Kondo physics within the projected 
Hilbert space of those four bands associated with the 
following operators: 

6 

Cnk = ^ ^AkaPAka and / ko . = f ka + ^2 W k/nk- (3) 

Act n—3 

Here, n = 1 — 6 labels the six p-orbital bands of i^MSB + 
i^soCo the state at k being denoted |v k ) = c^ k |0) with 
the coefficients ^ kcr . The localized /-electron operator 
is modified through Kondo hybridization with the higher 
energy bands, u£ a = (v£|Jf K |/ko-)/(M/ ~ <k)> ^ = 3-6, 
where is the energy of n-th band of Hmsb + Hsoc 
and |/kcr) = /kcr|0)- The effective Hamiltonian within the 
reduced Hilbert space of {|v k ), |v k ), |/k+), |/k-)}, where 

liU = fl\o), is 

n=l,2 a 

+ £ (r4i+^4 (4) 

n=l,2,cr 

The /-level dispersion and the effective Kondo cou- 
pling V£ a = (vg|iJ K |/ka) are given by 

j _ u , y IKI^K|/k.)| 2 
^ak + > 

VJT = J -f fa A - fcXk* - HtfJ) • (5) 

Some properties of the self-consistently obtained hy- 
bridization parameters of Kondo-coupled phases such as 
TJ = r£ = b z ao . = ¥ a ^ = b y ad . = are used to derive 
the above effective Hamiltonian. The energy dispersion 
obtained from H e R matches the full dispersion almost ex- 
actly. Characteristics of the Kondo-coupled bands can be 
discussed faithfully now in terms of the effective Hamil- 
tonian, i^eff- 

The effective Kondo coupling V£ a appearing in H e ^ 
has the explicit momentum dependence despite the fact 
that the original Hamiltonian only contained the on-site 
Kondo interaction. The symmetry of the Kondo pairing, 
self-consistently determined to be p-wave as illustrated 
in Fig. 4, also dictates that of the hybridized bands, 
as shown in Fig. 3. The symmetry propertes of V™ , 
in turn, follow from those of the wave functions ^ kcr as 
outlined in Table I. 
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wave function 
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l _ V x\i+ 5 _ V ?/k+ 5 V z\t+ 5 ^Xk- 5 V ?/k- 5 _ ^k- J 


( ky , k x ) 


{~^k+) ^Xk+ 5 ^k+ 5 —i V y\l- 5 ^Jk- 5 *^2k- } 



TABLE I: Relation between eigenstates of Hmsb + Hsoc at 
several symmetry-related momenta, in the first row 

refers to the amplitude of pA-orbital with spin a for the 
eigenstate at k = (k x ,k y ); a = +(— ) corresponds to spin 
up (down). Eigenstates at other momenta are related to this 
one as shown in the subsequent rows. 



(a) a=0 



(b) a=0.5 




FIG. 4: (Color online) Contour plots of |V£ + |, the absolute 
value of the effective Kondo hybridization, in the Brillouin 
zone, (a) a = (no SOC) and (b) a = 0.5. 



Transformation rules for the wave functions are de- 
duced by examining whether the Hamiltonians at a pair 
of symmetry-related momenta can transform into each 
other by some kind of unitary or anti-unitary opera- 
tors such as Ka z R x (tt), KI Ryi^), g z R z (tt) 
and e- i i 7T D 1 z /2 (3ir/2) ® R z (ir/2). We denote R x (0) for 

1/2 

the SO(3) rotation matrix, A = x,y,z, D z (0) for the 
SU(2) rotation matrix, and K for the complex conjuga- 
tion. From the Table and Eq. (5) one can show that 
y n Z h — —vVh^h 5 and therefore a unitary transfer- 
mation cr z <g> I between H G s(k) and i^ e ff(— k) exists, en- 
suring the equivalence of the energies between k and — k. 



and V n Z h are not 



On the other hand V n Z h , V^ a _ h 

related to V£ a k in any way according to the Table I. For 
example, V™J y = M-b^v%j2+h%v™ a -hjv™ a )/2 
does not even have the same magnitude as V£ ,<7 in gen- 
eral because the wave function i?™ ko . are complex- valued. 
Referring to Eq. (5), neither the reflection about the x- 
or the 7/-axis nor the 90° rotation in k-space is guaranteed 
to yield unitary-equivalent i? e ff- 

It is further possible to find an explicit expression of 
the effective Kondo coupling to prove its p-wave nature 
for models without SOC. Given the parameters such as 
V\ = 2, V2 = 0.65, 7 = 0.3 and the filling of conduc- 
tion electrons n c = 1.5, the band dispersions around 



the M point meet the conditions assumed in deriving 
the effective Hamiltonian i^ e ff- Here one can show ^ kcr 



near k = (7r,7r) 



is v. 



x,y\t+ 



-37s x>1/ /2(Vi + V 2 ) and 



J x,y\c- 



0, and vice versa for n = 2. Turning off SOC, 
all spin-mixing mean field parameters b z ad . vanish. We 
can conclude only two equivalent ground states with real- 
valued Kondo parameters h x = h y a or h% = —h y are 
possible, giving out 



3iJKh x j 
Wi + V 2 ) 



(sin k x ± sin k y ) (6) 



and V£ 



for h x = respectively. The ef- 
fective Kondo coupling indeed has the p-wave symmetry 
with a node along the k x = +k y or k x = — k y line. Figure 
4(b) shows the plot of |V^ + | for the case of h x = —h y 
where the node extends along k x = —k y . It is also seen 
that p-wave Kondo coupling is not predicated on SOC at 
all. 

Instead, the p-wave Kondo coupling is strictly a 
consequence of nonzero MSB, 7 / 0. When 7 = 0, 
Hamiltonian matrices of i^MSB + Hsoc at (k x ,k y ), 
(—k x ,k y ), (k x ,—k y ) and (—k X: —k y ) are equal to each 
other [13]. As a result, those momenta have identical 
eigenvectors and, according to Eq. (5), the same 
Kondo coupling V£ a . In addition the p-wave character 
arises only when both h x and h y a are nonzero. For 
phases I and II such condition is indeed achieved with 
\h x \ = I h v a I . The p-wave Kondo states are favored away 
from half-filling regime such as considered in the phase 
diagram of Fig. 2. At half filling n c = 3 without SOC, 
on the other hand, we often find the ground phase with 
only h z a non-vanishing. In this case, the effective Kondo 
coupling is in the form V£ a = —J¥.h%v^ a /2 w iJ^h%/2, 
which is s-wave as in the ordinary heavy fermion matter. 

Discussion - As the angstrom-scale deposition tech- 
nique of thin films develops into maturity, the role of sur- 
face states with broken mirror symmetry takes on greater 
practical significance. Such MSB field effect leads to in- 
teresting consequences when acting on a multi-orbital en- 
vironment such as the emergence of Rashba effect [6, 7], 
orbital-dependent spin transfer torque [14], and as argued 
in this work, of un-conventional Kondo pairing. Conven- 
tional one-band Kondo model with the phenomenolog- 
ical Rashba term could very well miss this feature due 
to the inadequate treatment of the multi-orbital char- 
acter. The asymmetric band structure we found in the 
Kondo-coupled phase should be readily detected in the 
spectroscopic study such as ARPES, or in the transport 
measurements through their directional responses. Al- 
though proposals of non-s-wave Kondo pairing are al- 
ready available [10-12], they tend to trace its origin to 
the non-local character of Kondo coupling between local- 
ized and itinerant moments [10], or the orbital symme- 
try of the /-orbital localized moments [12]. By contrast, 
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the p-wave Kondo hybridization emerging in our scenario 
is a generic consequence of the macroscopic breaking of 
mirror symmetry imposed by the geometry of the het- 
erostructure itself. Existing theory of heavy fermion het- 
erostructure [4] focused on the nature of unconventional 
superconductivity resulting from Rashba effect, but did 
not address the nature of Kondo pairing itself. 
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